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Abstract- A typical viscoplastic model will introduce up to three types of internal state variables in order 
to properly describe transient material behavior; they are: the back stress, the yield stress, and the drag 
strength. Different models employ different combinations of these internal variables — their selection and 
description of evolution being largely dependent on application and material selection. Under steady-state 
conditions, the internal variables cease to evolve and therefore become related to the external variables (stress 
and temperature) through simple functional relationships. A physically motivated hypothesis is presented 
that links the kinetic equation of viscoplasticity with that of creep under steady-state conditions. From 
this hypothesis one determines how the internal variables relate to one another at steady state, but most 
importantly, one obtains bounds on the magnitudes of stress and back stress, and on the yield stress and 
drag strength. 


EVOLUTION OF INELASTICITY 


The evolution of plastic strain, e?p in the classical theory of creep (ODQVIST [1974]) is given by 




(i) 


with the subscript implying steady state, and where ||e p [| denotes the magnitude of the plastic strain- 
rate, and ||S|| denotes the magnitude of the deviatoric stress, % The subscript signifying steady-state is 
not attached to Sij because stress is a controllable variable, whereas, creep-rate is a response variable and 
hence its attachment. This equation states that an increment in creep strain accumulates in the current 
direction of the deviatoric stress. A dot is placed over a variable to signify its time rate-of-change. 


The evolution of plastic strain in vis co plasticity (PRAGER [1949]) is defined as 


£P — i/ 2 ||e p ll 3 — 

,J ' 11 11 || S — B|| 


( 2 ) 


where Bij is the back stress which represents the anisotropic aspect of internal stress. This equation states 
that an increment in plastic strain accumulates in the direction of an effective stress, Sij— Bij. 


The norms, or magnitudes, pertaining to the deviatoric tensors of this paper are defined by 


||I|| = and \\J\\ = y/'/»JnJii, (3) 

where Iij is any deviatoric ‘strain-like 5 tensor, and Jij is any deviatoric ‘stress-like 5 tensor. These are the 
norms of von MISES [1913], where the coefficients under the radical signs scale the theory for shear. 

KINETICS OF EVOLUTION 


In the theory of creep, ||^ p ||^3 is described by a kinetic equation, i.e. an equation of state. ZENER k 
HOLLOMON [1944] determined that such a kinetic equation can, to a good approximation, be decomposed 
into a product of two functions; in particular, at steady state 


\\e p \ U, = mZss 


[Ml 

C 


> 0 , 


(4) 


where > 0 is a thermal function, Z > 0 is the Zener parameter, and C > 0 is a strength parameter that 
normalizes the stress. The Zener parameter is a temperature normalized measure of the plastic strain-rate. 
Square brackets, [*], are used throughout this paper to denote ‘function of’, and are therefore kept logically 
separate from parentheses, (•), which are used for mathematical groupings. 


The kinetics of viscoplasticity are also taken to be described by a ZENER k HOLLOMON [1944] decom- 
position of state, viz. (FREED et al . [1991]) 


= >o. 


(5) 


where the Macauley bracket, ((||S-B|| - Y)/D), has either a value of 0 whenever ||S — B|| < Y (defining the 
elastic domain), or a value of (||S — B|| - Y)/D whenever ||S — B|| > Y (defining the viscoplastic domain), 
with |)S— B|| = Y establishing the yield surface. The internal state variables — Bij , D and Y — are described 
by evolution equations that are functions of state. Remarkably, details of their evolution are not required in 
order for us to bound their values, with one exception, the back stress at steady state must be coaxial with 
the deviatoric stress so that the directions of plastic strain-rate defined by Eqns (1 & 2) are equivalent at 
steady state. 


BOUNDS ON INTERNAL STATE 


In order for a viscoplastic theory to reduce analytically to creep theory when at steady state (i.e. when 
B = 0, t> = 0 and Y = 0 for ||£ p || > 0) two conditions must be satisfied. First, the back stress at steady 
state must be coaxial with the deviatoric stress, as stated above. And second, it is necessary that the kinetics 
of viscoplasticity reduce analytically to the kinetics of creep unde^ steady-state conditions. To satisfy this 
second constraint, I hypothesize relationships between the steady-state and transient Zener parameters, and 
between the internal and external variables, when at steady state (FREED k WALKER [1990]). 


In support of the second constraint, given to secure equivalence between creep and viscoplasticity when 
at steady state, I shall suppose that (FREED k WALKER [1993]) 


Z — Z ss 


7 ||s-B||-r y 


( 6 ) 
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This relationship implies that the transient Zener parameter, Z, has the same functional form as the steady- 
state Zener parameter, Z 33 , but with a different argument. 

Furthermore, to also support this second constraint, I shall suppose that - 

I|B||— = / t ss [||S||] ||S|| , D a , = D 0 + 6\\S\\ and Y„ = (1-/)* M [||S||] ||S|| , (7) 

where i 38 > 0 and <5 > 0 are the steady-state fractions of applied stress that are associated with the 
internal stress (i.e. the back and yield stresses) and the drag strength, respectively, such that 0 < i < 1. 
The parameter / partitions the internal stress between isotropic and kinematic contributions, such that 
0 < / < 1. The fact that the drag strength is taken to be proportional to the applied stress at saturation 
is consistent with TAYLOR’S [1934] concept of a material’s innate strength to resist plastic flow, i.e. D is a 
strength parameter — not a stress parameter. We take the internal stress to be a nonlinear function of the 
applied stress at saturation because that is what the experimental data of ARGON & TAKEUCHI [1981] 
and CADEK [1987] suggest. 


Because the applied stress and the back stress must be coaxial at steady state, as discussed above, it 
follows that 

||S-B||„ = ||S||-||B|| m . (8) 

Therefore, upon equating arguments of the Zener parameters given in Eqns (4 & 5) in accordance with the 
hypothesis of Eqn. (6), while utilizing Eqns (7 & 8), one obtains the result 


C- Do -6 ||S|| C - Do + y/{C - A.) 2 - 4 6 C(||B||„ + Y„) 

1,5 ~ C 2 C ' W 

If one uses Eqn. (7) and writes ||B|| sg +y sg = t M [||S||] ||S||, then from Eqn. (9) one determines that ||B|| gg + 
Y aa = (C - Do - <5||S||)||S||/C. Because d(||B|| gg +y gg )/3||S|| = 0 establishes the maximum state of interned 
stress, one is lead to the result 

l|S||™ = ^. (10) 


Substituting this relationship back 

into Eqns (7 9) gives the additional upper bounds for: 

the back stress, 


urn f ( c ~ D °) 2 

| j-O ||ma* — / ^SC ’ 

(ii) 

the drag strength, 

Dmax — V2 , 



(12) 

and the yield stress, 




Y -ft f) (C~W 
r max -(i J) 4SC ■ 

(13) 


It is a remarkable fact that one can bound the stress and internal state variables without specifying the details 
of how these internal variables evolve. Similar bounds are given in FREED & WALKER [1993] for the case 
where the internal stress is composed of short- and long-range back stresses with no yield stress. 


Restricting i S3 to be real valued, and considering imin to be associated with the maximum attainable 
magnitude of internal stress, one finds that on approaching the limit of zero stress the ratio of interned stress 
to applied stress at steady state is at its maximum, i.e. 


v C ~ D » i 

lim C 33 = Lmax = ^ ^ 1 , 

l|S||-o c 


(14) 


3 



which is in reasonable agreement with the experimental observations of ARGON k TAKEUCHI [1981] and 
CADEK [1987]. Approaching the limit of maximum stress, this ratio attains its minimum, i.e. 


lim 

l|S|MlS|| 


= f'min — 


C-Do 

2C 


~ 1/2 


2 , 


(15) 


which is in reasonable agreement with the experimental observations of LOWE k MILLER [1983] and 
ARGON k BHATTACHARYA [1987]. 


To be physically meaningful, ||B|| > 0, D > 0 and Y > 0. Furthermore, their steady-state values ought 
to increase monotonically with increasing stress (Freed and Walker, 1990). This is verified easily for my 
hypothesis, Eqns (6, 7 k 8), as long as 0 < ||B|| ^ ||B|| moi , -Do ^ D ^ D ma x and 0 < Y < Y max • 
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